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In a simulation of SU(2) gauge theory we investigate, after maximal Abelian projection, the dual Maxwell 
equations for colour field and monopole current distributions around a static quark-antiquark pair QQ in vacuo. 
Within the dual superconductor picture we carry out a Ginzburg-Landau type analysis of the flux tube profile. 
As a result we can determine the coherence length of the GL wave function related to the monopole condensate, 
£ = .25(3) fm, to be compared to the penetration length, A = -15(2) fm (scaled with the string tension). 



1. INTRODUCTION 

The QCD vacuum can be viewed as a dual su- 
perconductor characterized by a monopole con- 
densate 0. When embedding a static QQ pair 
into the vacuum the latter expells by the dual 
Meissner effect the colour electric fiuxlines, thus 
giving rise to colour confinement. The core of 
the flux tube is just a normalconducting vortex 
which is stabilized by solenoidal magnetic super- 
currents, k, in the surrounding vacuum. 

These features are to hold in some effective 
Abelian subsector of the theory. In a series of pre- 
vious studies of SU(2) and SU(3) gauge theories 
it has indeed been confirmed that the important 
d.o.f.'s driving confinement are largely singled out 
by maximal Abelian gauge (MAG) projection [Q, 
as evidenced e.g. by the string tension 

In this note we will demonstrate that the va- 
lidity of this dual superconductivity picture can 
be extended to colour electric field and monopole 
current distributions, E and k, in rather de- 
tail. Refining previous studies along these lines 
|5|-0, we shall perform an unbiassed analysis of 
the Ginzburg-Landau (GL) wave function ||, /, 
which can be extracted in a sequence of steps from 
the profile of the flux tube in the transverse center 
plane (CP) to QQ, as seen on the lattice. / car- 
ries a second length scale (in addition to the dual 
photon mass, A), in form of the coherence length, 



£, of the condensate; it is the counterpart to the 
Abelian Higgs mass in the effective GL approach. 

2. EVIDENCE FOR DUAL THEORY 

It is crucial to start out by establishing the 
validity of the dual Maxwell theory. We create 
SU(2) field configurations in the usual way, go 
through MAG projection into the Abelian world 
Q and determine correlators between large Wil- 
son loops and plaquettes therein while k is 
defined according to the prescription of our lo- 
cal host |l0|. In Fig. [j] we verify that source and 
sink of E are strictly localized to the positions 
of Q and Q\ Fig. S illustrates that the data on 
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Figure 1. div E is strictly localized to the position 
of Q and Q, at separation r w 1.2 fm. 

CP transverse distributions (plotted versus radial 
distance in lattice units, x, out of core) provides 
impressive support of the dual Ampere Law. 
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Figure 2. Test of the dual Ampere Law, k = 
curlE. 

3. GINZBURG-LANDAU ANALYSIS 

Being nonlinear, the GL equations Q are not 
at all trivial to solve on a grid: 



m 3 +e 
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x dx V dx 



fix), 
(1) 



with 

m = (§-^) ■ < 2 ) 

Our approach is to determine / directly from 
the electric vector potential Ag (through integra- 
tion of E = curlA) and kg via the second GL 
equation (GLII). For this purpose we employ a 
parametrization for E z that respects the bound- 
ary conditions of / (BCf) which are rather in- 
volved, as they imply the value of E z (0) JTT] | . 

Given the functional form of E z the dual 
Ampere Law induces a parametrization for kg as 
well, which at the end leads us to a 4 param- 
eter ansatz to fit both distributions simultane- 
ously. For a discussion of the details we refer to 
ref. [jllj . The two interesting physics parameters 
resulting from the joint fit to E z and kg, namely 
the total colour electric flux $ (expected to have 
the value one) and the penetration length A come 
out to be $ = 1.08(2), A = 1.84(8), in consis- 
tency with the value from a deep London limit 



Figure 3. Fitting the azimuthal magnetic super- 
current and longitudinal colour electric field. 

analysis, A = 1.82(7) jLlJ. The simulations are 
performed on a 32 4 lattice at j3 = 2.5115, the 
lattice spacing being a « 0.081 fm ||. 

The quality of the fit is illustrated in Fig. ||, 
where we see nice agreement apart from some lat- 
tice artifacts in k in the small x regime, x < 2.2^|. 
As a result we would expect to be able to de- 
termine / from solving GLII outside this region. 
The result is shown in Fig. ^ where we plotted the 
CP transverse distribution of the wave function, 
together with the profile of the longitudinal field 
E z {x). Notice that we attain reasonable statisti- 
cal accuracy on /, up to about 4.2 lattice spac- 
ings out of the vortex core. Beyond this point, we 
are loosing overly in sensitivity for /. Neverthe- 
less one can observe / to approach asymptotia, 
/(oo) = 1, fairly rapidly. 

In order to finally extract the quantity of inter- 
est, £, it is useful to parametrize / with a func- 
tional form that again obeys BCf. To this end 
we try the usual ansatz f(x) — tanh(x/a), which 
is very successful, as evidenced in Fig. |[ with 
a = 3.33(5), but A = 1.62(2)0. Now that we have 
a good continuum form of f(x) from GLII at our 
disposal, we can plug it into GLI and solve the lat- 



2 We meet non-integer lattice spacings as we use off-axis 
separations! 

3 Note that a is not identical to the coherence length as 
sometimes assumed. 
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Figure 4. GL wavefunction and colour electric 
field in transverse profile. 

ter in terms of an effective function £(x). To the 
extent that this function proves to be constant we 
can claim to have arrived at a consistent deter- 
mination of the coherence length, £. In Fig. [| we 
exhibit the outcome: indeed £ shows only a rather 
mild (±10%) x dependency within our window of 
observation 2.2 < x < 4.2. Finally, translating 
the remaining variation into the systematic error, 
we can quote the value £ = 3.10l|g. 




4. DISCUSSION 

If one includes the uncertainties in the deter- 
mination of a and A, one obtains the error bands 
for £ as indicated in Fig. 0. A conservative esti- 
mate of the uncertainty in the determination of 
the penetration length, A = 1.84^24, blurs the es- 
timated transition point between type I and type 
II superconductivity, n = A/£ = 1/ \/2 into an er- 
ror interval, as indicated in Fig. ^ by the symbol 
"?" . We conclude that there is weak evidence for 
type I superconductivity. 

It would be worthwhile to confirm this tenta- 
tive conclusion by a scaling analysis and to study 
interaction among vortices by putting more than 
one of them into the vacuum. 
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Figure 5. Effective £ with error bands related to 
uncertainties in a and A. 



